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Abstract 

The consequences of recent experimental finding that hydrogen bonds of the anti- 
parallel /3-sheet in nonspecific binding site of serine proteases become significantly 
shorter and stronger synchronously with the catalytic act are examined. We investi- 
gate the effect of the transformation of an ordinary hydrogen bond into a low-barrier 
one on the crankshaft motion a peptide group in the anti-parallel /3-sheet. For this 
purpose we make use of a realistic model of the peptide chain with stringent mi- 
croscopically derived coupling interaction potential and effective on-site potential. 
The coupling interaction characterizing the peptide chain rigidity is found to be 
surprisingly weak and repulsive in character. The effective on-site potential is found 
to be a hard one, i.e., goes more steep than a harmonic one. At transformation of 
the ordinary hydrogen bond into the low-barrier one the frequency of crankshaft 
motion of the corresponding peptide group in the anti-parallel /3-sheet is roughly 
doubled. 
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1 Introduction 



The problem of localization and storage of energy by proteins at binding of 
ligands or protein-protein interactions is an old one in biophysics (for review 
see |1], |2] and refs. therein). This problem is especially urgent within the 
context of enzymatic catalysis. How does an enzyme store and utilize the 
energy released at substrate binding? The answer to this question is very 
poorly understood not only at quantitative but even at qualitative level. In 
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our opinion recent experimental data on serine proteases [3] are able to shed 
light on this problem. 

Indeed the high resolution X-ray diffraction studies (< 1.2 A) of serine protease 
intermediate structures revealed that "the strength of the hydrogen bonds 
between the enzyme and the substrate changed during catalysis. The well- 
conserved hydrogen bonds of antiparallel /3-sheet between the enzyme and the 
substrate become significantly shorter in the transition from a Michaelis com- 
plex analogue ... to an acyl-enzyme intermediate . . . presumably synchronously 
with the nucleophilic attack on the carbonyl carbon atom of the scissile pep- 
tide bond. This is interpreted as an active mechanism that utilizes the energy 
released from the stronger hydrogen bonds to overcome the energetic barrier 
of the nucleophilic attack by the hydroxyl group of the catalytic serine." [3j. It 
is worthy to note that these data are in coherence with those of (see also |5] 
and refs. therein) that the /3— strand corresponding to residues 214 — 217 shifts 
~ 0.8 A toward the inhibitor upon its binding by a— lytic proteases. Here one 
should strictly discriminate the established experimental fact of the hydrogen 
bonds shortening and strengthening in the antiparallel /3-sheet from the spec- 
ulation about the participation of this effect in the mechanism of the catalysis. 
The latter is an intriguing hypothesis that still needs to be confirmed. It will 
not be touched upon in the present paper. Our concern here is in modeling 
the consequences of the hydrogen bonds shortening and strengthening in the 
antiparallel /3-sheet for its dynamics. 

Serine proteases are the most extensively investigated enzymes (for review 
see [6], [7] and refs. therein). The greatest wealth of structural and kinetic 
data has been obtained for them. That is why these "working horses" of enzy- 
mology are the most suitable objects for investigating unresolved problems of 
enzymatic catalysis. The antiparallel /3-sheet mentioned above forms the non- 
specific binding site of serine proteases and is quite similar for all of them. The 
differences are in specific binding sites which are of no interest for us here. The 
antiparallel /3-sheet of the nonspecific binding site comprises three sequential 
peptide groups of enzyme backbone which form hydrogen bonds with corre- 
sponding peptide groups of the substrate (see Fig. 3 in [6] or a schematic picture 
Fig.24-9 in |H] or that Fig.11.6 in [H]). Despite the fact that two peptide chains 
belong to two different proteins the resulting structure has all features of the 
antiparallel /3-sheet and is usually called so without reserve. 

Nature of hydrogen bonds in proteins is a matter of controversy from the 
beginning of protein science and enzymology. In particular the possibility of 
the so-called low barrier hydrogen bonds in enzyme active sites is of utmost 
interest. The proposal about their role for enzymes was put forth by Gerlt, 
Clealand and Frey in 1993-1994 (for review see [10], [11] and refs. therein). 
Since then their participation in transition state stabilization is a matter of 
heated debates [T^j, [IS] - Especially that between His and Asp in the catalytic 
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triad of serine proteases remains to be the subject of numerous studies and 
non-stopping controversy with arguments pro [Il],[15], [16], [T7], [18], [T9] . 
[20] . [21], [22] and contra [5], [13], [23]. We do not touch upon the catalytic 
triad in the present paper and focus ourselves exclusively on the processes 
in the nonspecific binding site. The most salient feature of low barrier hy- 
drogen bonds as compared to ordinary hydrogen bonds is their strength and 
short distance between hydrogen donor and acceptor atoms (dDA). We do not 
involve ourselves in very intricate classification debates about blurred bound- 
aries between weak, strong, short strong, short ionic, low barrier and very 
strong hydrogen bonds (see, e.g., [21] )• For convenience we freely divide the 
hydrogen bonds to strong ones (SHB) and weak ones (WHB). We set that 
SHB have energies > 10 kcal/mol and dDA < 2.5 A while WHB have energies 
2-^5 kcal/mol and dDA > 2.8 A. Thus the aim of the present paper is to 
investigate the transformation of WHB into SHB in the antiparallel /3-sheet of 
the nonspecific binding site of serine proteases. To be more precise we study 
the effect of this transformation on the backbone dynamics in the antiparallel 
/3-sheet. However it should be stressed that following [23] we believe that there 
are no well defined boundaries between different types of hydrogen bonds. In 
fact we investigate the parametric dependence of the backbone dynamics on 
the hydrogen bond energy from the range 2 ^ 30 kcal/mol. 

Last years are marked by noticeable progress in revealing subtleties of pro- 
tein dynamics gained by infra-red (IR) spectroscopy [25], [26], [2Z], [21], [29] . 
[30] . normal mode analysis [31], [32], [33] and molecular-dynamics simulations 
(see the recent reviews [33], [32] and refs. therein). There are several analyti- 
cal models treating mechanics of peptide chain backbone [36], [3Z], [33], [39] . 
The most interesting type of motion for our context is the so the so called 
"crankshaft-like" one. It means the rocking of the rigid plane of the peptide 
group due to degrees of freedom of torsional (dihedral) angles tpi and ilJi-i (see 
Fig.l). It was proposed on theoretical grounds from normal- mode analysis [30] 
and is supported by NMR experiments and molecular dynamics simulations of 
protein backbone [31], [32], [33], [33], [3S], [3S], [3Z], [33], [3S], [SO], [51], [S2], 
[53], [53]. It is comprehended now as a dominant type of motion for the latter 
that "involves only a localized oscillation of the plane of the peptide group. 
This motion results in a strong anticorrelated motion of the $ angle of the 
z— th residue ($i) and the ^ angle of the residue i — 1 (^E'i-i) on the 0.1 ps 
time scale." [33] ■ Thus the essence of this motion is the so called anticorrelated 
motion of the torsional angles and ipi-i manifested itself in the requirement 
(see Fig. 2) 

0/2 = A^, = A^i_i 

In this case the plane of the peptide group rocks as a whole around some axis 
a that goes through the center of masses of the peptide group parallel to the 
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bonds C^~^ — and N'^ — C^, (see Fig.l). The moment of inertia of the peptide 
group relative to the axis a can be easily calculated to be / ~ 7.34:-10~^^g-cm^ 
(see Appendix). Molecular dynamics simulations and NMR experimental data 
suggest that the character of the correlation function for the crankshaft motion 
is decaying oscillations [13] , [51] but provide characteristics for them in a very 
wide range from the subpicosecond and picosecond time scale [51] to slower 
motions on a much larger time scale from tens of picoseconds to 100 ps and 
more [53J . This is presumably a manifestation of the " Russian doll" structure 
of the conformational potential for the crankshaft motion when a group of local 
minima forms a smooth local minimum and so on. In fact knowing the actual 
values of the frequency for the oscillations and the characteristic time of their 
decay for the functionally important crankshaft motion is not indispensable for 
the purposes of the present paper. However in our opinion it is reasonable to 
assume that the frequency of oscillations of the plane of the peptide group as 
a whole for such motion should be at least order of magnitude less than those 
for high frequency in - plane motions such as, e.g., Amide-I (~ 1600 cm'^). 
The choice of the frequency in the uq ~ 100 cm~^ ~ 10^^ range enables 
us to match it with the amplitudes of rocking of order of several degrees (see 
(fT^ below) in accordance with experimental data [5U|, [52] . 



To investigate the effect of the transformation of WHB into SHB in the an- 
tiparallel /3-sheet on the characteristics of the backbone crankshaft motion we 
use the model developed earlier [55j, [56]- The latter is a microscopically strin- 
gent model of polypeptide backbone dynamics. It deals with realistic effective 
on-site potentials and coupling interaction potentials but at the same time 
is simple and tractable. The model yields for the correlation function of the 
crankshaft motion the exponentially decaying oscillations [SB] in qualitative 
agreement with the results of molecular dynamics simulations [M] ■ The model 
enables us to elucidate the role of hydrogen bonds in the backbone dynamics. 



The paper is organized as follows. In Sec. 2 the Hamiltonian of the model is 
constructed. In Sec. 3 the coupling interaction of two adjacent peptide groups 
is considered. In Sec. 4 the local on-site potential for the peptide group is con- 
sidered. In Sec. 5 the contribution of hydrogen bonds is separately considered 
because they are the central point of the present paper. In Sec. 6 the dimen- 
sionless characteristics of the effective on-site potential and of the coupling 
interaction potential are discussed. In Sec. 7 the equation of motion is derived 
and the correlation function for crankshaft motion is obtained. In Sec. 8 the 
results are discussed and the conclusions are summarized. In Appendix some 
technical details are presented. 
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2 The Hamiltonian of the model 

A schematic picture of a polypeptide chain with all designations used further 
is presented in Fig.l. The mutual orientation of two adjacent peptide groups 
is characterized by the torsional (dihedral) angles and ipi. The angle ipi 
characterizes the rotation round the bond Cf — Ni and that ipi characterizes 
the rotation round the bond — Q. Further we consider the equilibrium 
dynamics in the antiparallel /3-sheet for which the equilibrium values of the 
torsional angles are the same for all peptide groups (they are listed at the end 
of the next Sec). We consider small deviations (pi(t) and 'i/'i(^) of the torsional 
angles from their equilibrium values {ipi = + (pi{t) and ipi = ip^ + ipiit)) 
with |v2i(t)| < 20° and |V'j(^)| < 20°. At such amplitudes of the deviations the 
hydrogen bonds confining the peptide group in the antiparallel /3-sheet are not 
broken [8]. Finally we consider a peculiar type of motion ilji_i{t) = fiit) that is 
namely the crankshaft one. The latter means that the peptide group rotates as 
a whole respective some effective axis a passing through the center of the C— N 
bond parallel to the bonds Cf — Q and A^j+i — C°|_i that are assumed to be 
approximately parallel {ZC^-C'i-Ni+i = 113° and ZC[- Ni+i-C^^^ = 123° 
[57]). This type of motion is stipulated by the fact that the peptide group is 
a planar rigid structure [8], [S7]. For the sake of uniformity of designations we 
further denote Xi = 2ipi(t) = 2ipi_i{t) (see Fig. 2) and thus 



The moment of inertia of the peptide group relative to the axis a can be easily 
calculated and is J ~ 7.34 * IQ^^^g * cm^ (see Appendix). The Hamiltonian of 
the polypeptide chain in our model with nearest neighbor interactions is 



Here Uioc{xi) includes interactions defining the local potential of the peptide 
group (namely hydrogen bonds, the so-called torsional potentials and the van 
der Waals interaction of covalently non-bonded atoms of the peptide group 
with the atoms of adjacent side chains) defining its separate motion while 
U{xi; Xi^i) includes the coupling interactions (namely the van der Waals in- 
teraction of covalently non-bonded atoms of adjacent peptide groups i and i-|-l 
and their electrostatic interaction) intermixing the motions of these groups and 
leading to the rigidity of the peptide chain. It should be stressed that the lat- 
ter potential also contributes into the separate motion of the peptide groups. 
We can define the effective on-site potential for such motion as 



€-1 + ^^/'^ 



(1) 




(2) 



Veffixi) = Uiocixi) + U{xi] 0) + t/(0;xi) 



(3) 
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and the coupling interaction potential as 



Thus the equation of motion is 

^cPxi _ dUiocjxi) _ dU{xi^i]x,i) _ dU{xi] Xj+i) 
diij doc 2, doc ^ doc ^ 



(5) 



In the following two sections we consider in details the functions Uioc{xi) and 
Xj+i). At doing it we will freely pass back and forth between the variables 
according to the rule ([1]) which can be rewritten as 

Av3j=Xi/2; A?/'i = Xi+i/2 (6) 



3 Coupling interaction defining the rigidity of a polypeptide chain 



Both interactions in (j4]) intermixing the motions of the adjacent peptide groups 
and contributing to the coupling interaction are described by the central po- 
tentials W^^{Rmn{,'^i] ■ipi)) and W^^{Rjnn{y^i; ■ipi)) between the atom A„ of the 
i-th peptide group and the atom of the i + 1-th one with R^n being the 
distance between the atoms An and A^- The electrostatic potential is 

W^^'. (Rmn) = ^ (7) 



where and g„ are partial charges on the atoms A^ and An respectively 
[q{N) = -0.28e; q{H) = 0.28e; g(0) = -0.39e; q{C) = 0.39e where e = 
4.8* 10^^*^ CGS [57]) and e is the dielectric constant which for protein interior 
should be better conceived as some adjustable parameter 2 10 with 
3.5 being commonly accepted value close to high frequency permeability of 
peptides) |57j, [8]. In some studies e is supposed to be solvent dependent and 
chosen, e.g., 4 for CCh, 6^7 for CHCh and 10 for H2O [58], |59]. For the 
van der Waals potential one can choose any of the numerous forms suggested 
in the literature, e.g., the Lennard- Jones one (6-12) or the Buchinghem one 
(6-exp). In this chapter we use for numerical estimates the former one 

A R 

W [Kmn) = + l«J 

with the set of the well known parameters of Scott and Scheraga (other sets 
were also verified and found to give similar results). 
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The distance Rmn as a function of the angles (pi and ijji is 

Rmn{y^i,1pi) = {p'L + rl + 

2pmrn[sin 6'(cos7m sin a„ cos — sin7m cosa„ cosipi) — 

sin 7^ sin (cos ^/'j cos (fi cos 6 — simjji sin (fi) — cos 7^ cos q;„ cos 6*]} ^''^(9) 



The interaction potential for both van der Waals and electrostatic interactions 
has the form 



{vdw \ f vdw \ 

. i{xf,Xi+,)=J2Wl'' f{Rmn{Xi;Xi+,)) (10) 



where the summation in n is over atoms in the z-th peptide group and that in 
m is over atoms in the i + 1-th one. 

In what follows we work with full realistic coupling interaction potential de- 
scribed above. However for making use of suggestive analogies with the known 
literature results obtained on toy models we calculate the so called coupling 
constant which is a key characteristic of the truncated coupling interaction 
potential. Expanding the potential we obtain that to the leading order in the 
terms <^ 1 and Xj+i ^ 1 the rigidity of the peptide chain is determined by 
the term of the potential which intermixes the motion of the adjacent peptide 
groups 

U'^'^{xf,x,+,)^{-K)x,Xi+, (11) 



where the coupling constant is 
-K = -{K"'^'^ + K^^) = 

d^U"^^ (x,;a:,+i) ^ d'U^' (x,;x,+i) ^^^^ 

dXidXi+i |a;,=0;xi+i=0 dXidXi+i |xi=0;Xi+i=0 



The results of calculations of the contributions into the coupling constant for 
the antiparallel /3-sheet are as follows: v?° = -139°; = 135° [57]:-/^^'^"' ^ 
-5.0 ■ 10^1^ erg; -K^^ ^ 2.64 ■ 10"^^ erg; -K ^ 2.14 ■ lO^^^ erg aX e = 3.5. 
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4 Local on-site potential of the peptide group 

The local potential Uioc{xi) is composed as follows 

Uloc{Xi) = Usc{Xi) + Uhb{Xi) + Utors{Xi) (13) 

Here Usc{xi) is the energy of van der Waals interactions of the atoms of the 
i-th. peptide group with those of the side chains Ri and Ri-i and also with the 
atoms Hl^ and Uhb{xi) = u'f^ {xi) + U^^ {xi) is the energy of two hydrogen 

bonds of the 2-th peptide group and Utors{xi) = Uforsi.'^i) + Uf^rsi'^i-i) is the 
energy of the torsional potentials for the rotation of the angles ipi and ipi-i- 
For the latter we take the usual form 

Uforsi.'^'i) = E^{1 + COs3v3i) 

[/t.(V^^) = ^^(l + cos3V^,) (14) 

with E^p ^1 kcal/mol and ^1 kcal/mol. In all further simulations we choose 
the side chain to be Ala {Ri is Cp{H)^) and initially find the value of the angle 
Xi (see Fig.l) to minimize the energy. 



5 Hydrogen bond potential 

In f|T3l) uj£{xi) (j=l,2) are the contributions of two hydrogen bonds which 
a peptide group in the anti-parallel /3-sheet forms with its partners. We start 
from a model Morse potential usually employed to describe the hydrogen bond 
energy dependence on its current length [58], [59] 

Uhbix^) =D{1- exp [-n{r{x,) - ro)])' - D (15) 

Here n ~ 3 and tq ^ 1.8 A are adjustable parameters for distance de- 
pendence and D is an adjustable parameter of the hydrogen bond energy. For 
WHB D 4 6 kcal/mol at e = 3.5 with D = 5 kcal/mol being a conven- 
tional value [8] (in [59] D is assumed to be a decreasing value at increasing e 
with D = 0.5 kcal/mol at e = 10). Finally r{xi) is the current length of the 
hydrogen bond. 

We denote a^^^ the distance from the axis a (see Fig. 2) to the carbonyl oxygen 
atom (j = 1 so that a^^^ 1.5 A (see Appendix)) or that to the hydrogen 
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atom (j = 2 so that a*^^^ ^ 2 A (see Appendix)). Then the current length is 
related to the angle by straightforward trigonometry 

(rO-) (xP))' = (/(^■))' + 4 (a(^))' (sin(xpV4))' [a^^'^ + l^'^ cos 9^^^ 



This expression takes into account that the range for the angular displacement 
Xi (defined in ([T])) is in because it is twice of the torsional angles ipi and 
(see ([1])). However we also want to take into account that the current 
length decreases with the increase of the hydrogen bond energy. We attain it 
by introducing a linear dependence factor (1 — qD) in the expression for the 
current length where q is an empirical parameter 



0) ^U) 



qD 



X 



y + 4 {aU)f (sin (xpV4))' [a^''^ + l^^^ cos Oi^)] 



Thus for the f/ii^(x?'^) (j=l,2) we obtain 



C/(J(a;p)) = -i^(^■)+Z}W<!l- 



— n — r. 



{l-qD 



X 



y (/O'))' + 4 (aO'))' (sin (x?V4))' [a^^'^ + l^^^ cos ^0')] 



n ^ 2 



(16) 



Here n = 3 A^^ and ro = 1.8 A are old adjustable parameters while q and 
l^^^ are two new ones. The latter can be evaluated from the following data 
obtained in [60]. A 0.5 A lengthening of the low-barrier hydrogen bond results 
in a weakening of that bond by over 6 kcal/mol. A 1.0 A lengthening of 
the hydrogen bond results in an approximately 12 kcal/mol decrease in the 
calculated strength of the corresponding hydrogen bond. Taking the energy of 
the low-barrier hydrogen bond to be D = 15 kcal/mol we obtain from these 
requirements and (|T6|) that q ~ 0.01 mol/kcal and l^^'^ ~ 2.0 A. The angles 
6**^^^ and 9^'^'^ can be evaluated from the known data on the bond lengths and 
angles in the peptide group ^7] as 9^^^ ^ 9^'^^ ^ 30°. 



6 Effective on-site potential and coupling interaction potential 



The adjustable parameters are chosen to satisfy the following set of require- 
ments: 1. The effective on-site potential Veff{x,i) (see [3]) must have a minimum 
at equilibrium {xi = 0) because the antiparallel /9-sheet is a steady stable 
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structure. 2. For the anti-parallel /S-sheet the hydrogen bond are known (see, 
e.g., [9] or [57]) to be practically parallel to the corresponding covalent bonds 
C = O and N—H (see Fig. 2). 3. For the anti-parallel /3-sheet the spectroscopic 
frequency 



K'}/(x^ 0)// 

27rc 



(1/A)sp = ^LJ1^____ (17) 



are ~ 100 cm ^. Here A is a wavelength, c is the light speed and the dash 
denotes a derivative in Xj. 

We take for the anti-parallel /3-sheet ^9° = —139°; ip'^ = 135° [57j. One hydro- 
gen bond for the peptide group the anti-parallel /3-sheet (namely that with 
j = 1) always remains to be WHB {D^^'' = 5 kcal/mol). For the second hy- 
drogen bond we vary the value D'-^^ for the energy within the range 2 30 
kcal/mol. 

In what follows we use dimensionless variables. We define the frequency 



and measure time in the units of ^ so that the dimensionless time is 

T = tu (19) 

Also we measure energy in the units of lu'^ and denote the dimensionless 
coupling constant p 



7 Equation of motion and correlation function for crankshaft mo- 
tion 



In this Sec. we denote for simplicity Xi{t) = (j){t) and rewrite ([T]) in the form 

0/2 = Aif, = AV^,_i (21) 

It is natural to describe the crankshaft motion by a Langevin equation. Such 
equations are frequently used in protein dynamics [6l], [62], [63], [M], [65] . 
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[66]. The equation of the crankshaft motion for the rigid plane of the peptide 
group in the above conditions is 



d'^(j){t) d(j){t) 



{22) 



where 7 is the friction coefficient and ^(t) is the random torque with zero 
mean < ^(t) = Q > and correlation function 



< mm >= 2kBT^5{t) 



(23) 



ks is the Boltzman constant, T is the temperature. 



We consider the hydrodynamic friction. This "macroscopic" notion is known to 
work surprisingly well at the molecular level (see [67J for thorough discussion). 
We model the peptide group by an oblate ellipsoid with half-axes a, h and c 
(where a ~ c and a, c >> h that reflects the fiat character of the peptide 
group) Fig.8. For the rotation of the ellipsoid around x-axis (that is in our 
case actually the axis a for the peptide group introduced above) the friction 
coefficient is given by a formula 



{l)x = 8-717] abc 



12r 



+ 



{be 



1 1/4 



,1/2 



(1 + ^ 



4r 



/2(6c)l/4 



(24) 



where r = 2(a6c)^/^ and r] is the viscosity of the ellipsoid environment. The 
required behavior is obtained if we have the condition of the underdamped 
motion 
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« 1 



(25) 



At such requirement we obtain (see e.g. [69], [63] ) 



a{t) =< 0(O)0(t) 



■exp 



l\t 
21 



cos (uJot) 



(26) 



We denote 



1^ 



7 



V2IUJ0 



(27) 
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Then the Fourier spectrum for the correlation function a{t) is 



aiu] 



1 + /i^ + {co/uoY 



^ + 2^2 1 + {u/u^f + 1 - (w/wo) 



(28) 



At our requirement (l25l) it has a sharp resonance peak at the frequency ujq. 
For the mean squared amphtude (msa) we have 



msa = ^ „ (29) 



The latter means that the amplitude of the crankshaft motion at room tem- 
perature and, e.g., typical value loq ~ lO^^s"^ is 

<Pma. = ^1^^0.1 ^6° (30) 



that is Aifi = Aipi-i ~ 3°. This value sheds light on the origin of essentially 
vibrational character of the peptide group motion manifested itself in fl26|) and 
fl28l) . At such angles of rotational deviation of the peptide group from its mean 
averaged position the linear displacements of the atoms are ~ c-cpmax ~ 0.1 A 
that is much less than both the size of a solvent molecule and interatomic 
distances to neighbor fragments of protein structure. That is why the envi- 
ronment exerts rather weak friction for such type of the peptide group motion 
that is reflected in the requirement (!25l) . Thus we conclude that the thermally 
equilibrium crankshaft motion of the peptide group is of essentially vibrational 
character even in such condensed medium as protein interior. Our model yields 
for the correlation function of the crankshaft motion the exponentially decay- 
ing oscillations f l26|) . This type of behavior is in qualitative agreement with 
the results of molecular dynamics simulations f5^. 



8 Results and discussion 



The effective on-site potential for the peptide group in the anti-parallel /3-sheet 
for the case of one weak hydrogen bond [D^^^* = 5 kcal/mol) and one strong 
hydrogen bond {D^'^^ = 30 kcal/mol) is presented in Fig. 4. In Fig. 5 the depen- 
dence of the spectroscopic frequency of crankshaft motion (1/A)sp (see ffT7|) ) for 
the peptide group in the antiparallel /3-sheet on the energy of one of its hydro- 
gen bonds (namely D^'^^) at different values of the torsional potentials for the 
rotation of the angles ipi and ipi^i (see (fT^ ) is depicted. This Fig. shows that 
at transformation WHB — )■ SHB of WHB with energy 5 kcal/mol into SHE 
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with energy say 15-i-20 kcal/mol the frequency of crankshaft motion becomes 
roughly twice of what it was ((l/A)^^ ~ 40 cm~^ — )■ (1/A)sp ~ 80 cm~^). In 
Fig. 6 the dependence of the dimensionless couphng constant p (see (|2Ui) ) for 
the peptide group in the antiparallel /3-sheet on the energy of one of its hydro- 
gen bonds (namely D^'^^) at different values of the torsional potentials for the 
rotation of the angles ipi and ipi-i (see f|T^ ) is depicted. The dimensionless 
coupling constant decreases in absolute value with the increase 

ofD(2) 

because 

it is inversely proportional to the u'^ (see f l20|) ) while the latter increases with 
the increase of D^'^^ (see Fig. 6). The dimensional value of the coupling con- 
stant K remains unchanged with variation of D^'^'> because does not depend 
on the hydrogen bond potential (see (fT2!) ). In Fig. 7 the dependence of the cou- 
pling constant p (see (l20l) ) for the peptide group in the antiparallel /3-sheet 
on the dielectric constant e at different values of the energy of one of its hy- 
drogen bonds (namely D^"^^) is depicted. This Fig. shows that the dependence 
of the coupling constant p (see (12U1) ) for the peptide group in the antiparallel 
/3-sheet on the dielectric constant e is rather strong. On the other hand the 
dependence of the frequency of crankshaft motion for the peptide group in the 
antiparallel /3-sheet on the dielectric constant e is negligibly small. This can 
be explained as follows. The frequency of crankshaft motion is determined by 
the the effective on-site potential (see [3] and [121) that does not include electro- 
static interactions while the coupling constant is determined by the coupling 
interaction potential (see [12] and [20!) that strongly depends on electrostatic 
interactions and thus very sensitive to the value of e. 

Our results of the investigation of the effective on-site potential and the cou- 
pling interaction potential for the anti-parallel /3-sheet testify the following: 
1. For the anti-parallel /3-sheet the coupling interaction is surprisingly weak 
I p |<< 1 (though the same is true for other types of protein secondary struc- 
ture [SS]). 2. The effective on-site potential is hard (goes more steep than a 
harmonic one). 3. The coupling interaction is repulsive p < at reasonable 
values of the dielectric constant e < 15 (the latter means that a non-zero value 
of a peptide group displacement tends to increase the values of the neighbor- 
ing peptide groups displacements with the opposite sign). The weakness of 
the coupling interaction | p |<< 1 is known to be indispensable for the possi- 
bility of creation and existence of a specific excitation (the so-called discrete 
breather) in a chain of non-linear oscillators. Thus we conclude that in the 
anti-parallel /3-sheet a discrete breather can be excited. 

The fact that dimensionless values of the coupling constant p are very small, 
i.e., I p |<< 1 (see Fig.6 and Fig. 7) means that there is little interdependence 
in the rocking (i.e.. Crankshaft motion Aipi = Aipi^i) of the planes of the 
adjacent peptide groups in the antiparallel /3-sheet. This result is in agreement 
with the conclusion of the authors of \5^ that "only a slight correlation is found 
between the motions of the two backbone dihedral angles of the same residue" 
(i.e., Aipi practically does not correlate with Aip^). 
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The most intriguing question is how the the transformation of WHB into SHB 
in the antiparallel /3-sheet of the nonspecific binding site of serine proteases 
can be related to the mechanism of enzyme action? The hypothesis that such 
relationship can take place was suggested in [5] (see Introduction). In our 
opinion it is very interesting and requires discussion and further development. 
The indispensable first step is the storage of the energy released at such trans- 
formation. The second step is the transmission of this energy to the scissile 
bond in the substrate molecule. There are several possible mechanisms of such 
storage in the antiparallel /3-sheet. One of them makes use of the phenomenon 
of the so-called discrete breather. This possibility was discussed in [HS], [TP] . 
The energy of the discrete breather can be transmitted to the scissile bond 
in the substrate molecule by the electric field fluctuations |70]. Another inter- 
esting possibility is the excitation at the WHB — )■ SHB transformation of the 
self-trapped Amide-I vibration in the C = O stretching mode of the nearest to 
the active center peptide group in the antiparallel /3-sheet of the nonspecific 
binding site of serine proteases. Indeed the Amide-I quantum energy is known 
to be 0.205 eV (see, e.g., [57]), i.e., ~ 4.7 kcal/mol. Thus the transformation of 

a WHB N — H O = C with energy 5 kcal / mol into a low-barrier hydrogen 

bond N — H •••O = C with energy, e.g, 15 kcal/mol releases enough energy 
to excite Amide-I of the C = O bond. Of course one should keep in mind 
that there can be complications due to the fact that the bond between C and 
O atoms in a peptide group is only partially double one (see, e.g., [57] or [9]). 
However these complications do not seem to be crucial. If the C = O vibration 
is excited then the energy can be transmitted to the scissile bond Ccarbonyi — N 
in the peptide group of the substrate molecule by the nonradiative resonant 
interaction of the dipole moment d of the stretching mode with that of the 
Ccarbonyi = O boud adjaccut to Ccarbonyi — N . At the Oversimplifying approx- 
imation that both dipoles have equal tilt relative the axis z connecting them 
the exchange energy (vibrational coupling) of the dipole-dipole interaction J 
can be evaluated as |7T], [72] 

J = -^f3cos20-l) 

Here R is the distance between the C = O bonds and Q is the angle between 
the dipole moment J and the axis z. The value of J for, e.g., peptide groups 
in a— helix was estimated from experimental data of infrared spectra and it 
was found that J = 1.55 ■ lO'^^ Joul [73], i.e., J ^ 2.22 ■ lO'^ kcal/mol. Thus 
the Hamiltonian for two interacting C = O bonds can be written as 

2 2 

Hex = ^0 X/ ^n^n ~ X/ i^^n^n+l + O'nO'n+l) 
n=l n=\ 

where ujq is the frequency of Amide-I (~ 1600 cm~^) and a^(a„) is the creation 
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(annihilation) operator for an Amide-I exciton in the site n. The storage of 
vibrational energy in the C = O stretching mode is widely used within the 
framework of the Davydov model for the bio-energy transport in a polypep- 
tide chain [7T], [72]. In this model the C = O vibration is excited by the 
energy released in the adenosine triphosphate (ATP) hydrolysis that yields 
about 0.43 eV. Hamiltonian of the above type is still a matter intensive inves- 
tigations within the context of energy transfer in proteins (see, e.g., [71] and 
refs. therein). Regretfully we failed to find experimental or theoretical investi- 
gations of the key point for our hypothetical scenario namely that of Amide-I 
excitation at transformation of the ordinary hydrogen bond into low-barrier 
ones. That is why we restrict ourselves from further discussing this speculation 
and leave the development of posed above possibility for future work. 

The conclusions of the paper are summarized as follows. The dynamics of a 
peptide chain in the anti-parallel /3-sheet is considered within a realistic model 
with stringent microscopically derived coupling interaction potential and ef- 
fective on-site potential. The coupling interaction is found to be surprisingly 
weak and repulsive in character. The effective on-site potential is found to be 
a hard one. At transformation of ordinary hydrogen bond into low-barrier one 
the frequency of crankshaft motion of the corresponding peptide group in the 
anti-parallel /3-sheet is roughly doubled. 

Acknowledgements. The author is grateful to Prof. Yu.F. Zuev and Dr. B.Z. 
Idiyatullin for helpful discussions. The work was supported by the grant from 
RFBR and the programme "Molecular and Cellular Biology" of RAS. 



9 Appendix A 

Here we give a brief estimate of the moment of inertia of the peptide group 
for rotation around the axis a that goes through the center of masses of the 
peptide group parallel to the bonds C^~^ — C and A^* — (see Fig.2). The 
peptide group is a rigid plane structure so that all atoms O, C, N and H 
lie in a plane. That is why the moment of inertia is simply the sum of their 
masses multiplied by the square of their distances to the axis cr. The masses 
in atomic units (1 a.u.m. = 1.7 ■ 10"^^^') are mo = 16, mc = 12, tun = 14 
and tuh = 1. That is why the center of masses is shifted a little to the atoms 
C and O relative the middle of the bond C — N. The lengths of the bonds 
are Iqc = 1-24 A, Iqn = 1.32 A and Inh = 1-0 A. That is why the distances 
to the axis a are approximately tq ~ 1.5 A, rc ~ 0.3 A and ^ 0.7 A and 
rjj ^ 2 A. Substituting these values into the formula / = Y^rriirf we obtain 
/ 7.6 ■ IQ'^'^g ■ cm? . Thus our rough estimate corroborate the precise value 
/ ^ 7.34 • IQ^^^g ■ cm^ from [75]. 
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Fig. 1. A schematic representation of the peptide chain with all designations neces- 
sary for the text. The covalent bonds are shown as bold lines. Dotted lines denote 
hydrogen bonds. Dashed lines denote auxiliary axes. 
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2. Schematic picture of the peptide group in the antiparallel /3-sheet . 
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Fig. 3. A look on the peptide group from the axis of rotation a explaining 
definition of the angle (defined in ([3]) ) . 
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Fig. 4. The normalized effective on-site potential Ve//(x) (see ([3])) for the peptide 
group in the antiparallel /3-sheet in case of one weak hydrogen bond (D^^^ = 5 
kcal/mol) and one strong hydrogen bond {D^"^^ = 30 kcal/mol). The value of the 
dielectric constant is e = 3.5. Here I is the moment of inertia of the peptide group 
relative the axis a and uj is the frequency. The range for the angular displacement 
X (defined in ([T])) is Att because it is twice of the torsional angles ip and ^ (see ([T])). 
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Fig. 5. The dependence of the spectroscopic frequency of crankshaft motion (1/A)sp 
(see pT|) ) for the peptide group in the antiparallel /3-sheet on the energy of one of its 
hydrogen bonds (namely D^'^^) at different values of the torsional potentials for the 
rotation of the angles ipi and ipt-i (see dHD): =0.5 kcal/mol, =0.5 kcal/mol; 
Ey, =1. kcal/mol, E^ =1. kcal/mol; E^p =1.5 kcal/mol, E^ =1.5 kcal/mol; E^ =2. 
kcal/mol, E^ =2. kcal/mol. The value of the dielectric constant is e = 3.5. 
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Fig. 6. The dependence of the couphng constant p (see (120)) ) for the peptide group 
in the antiparallel /3-sheet on the energy of one of its hydrogen bonds (namely D^'^^) 
at different values of the torsional potentials for the rotation of the angles ipi and 
Tpi-i (see p4p): E^p =0.5 kcal/mol, =0.5 kcal/mol; E^ =1. kcal/mol, E^ =1. 
kcal/mol; E^ =1.5 kcal/mol, E^ =1.5 kcal/mol; E^ =2. kcal/mol, E^ =2. kcal/mol. 
The value of the dielectric constant is e = 3.5. 



26 



0.000 



-0.005 - 



▼ D"'=30keal/mol 



A D'-'=20kcal/mol 



D'^' = 10kcal/mol 



D'^'=5 kcal/mol 



-0.010 



-0.015 - 



10 



12 



14 



Fig. 7. The dependence of the couphng constant p (see (f20]l ) for the peptide group in 
the antiparallel /3-sheet on the dielectric constant e at different values of the energy 
of one of its hydrogen bonds (namely D^"^^): D^'^^ =5 kcal/mol; D^^^ =10 kcal/mol; 
Z)(2) =20 kcal/mol; D^'^^ =30 kcal/mol. The values of the torsional potentials for the 
rotation of the angles (pi and tpi-i (see (jl4[) ) are =1. kcal/mol, =1. kcal/mol. 
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Fig. 8. Model of the peptide group by oblate ellipsoid. 



28 



